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^ Abstract 

We investigate and classify Fermi surface behavior for a set of fermionic modes in a family 
of backgrounds holographically dual to A/" = 4 Super- Yang-Mills theory at zero temperature with 
two distinct chemical potentials. We numerically solve fluctuation equations for every spin-1/2 
field in five-dimensional maximally supersymmetric gauged supergravity not mixing with gravitini. 
^ Different modes manifest two, one or zero Fermi surface singularities, all associated to non-Fermi 

liquids, and we calculate dispersion relations and widths of excitations. We study two limits where 
the zero-temperature entropy vanishes. In one limit, a Fermi surface approaches a marginal Fermi 
liquid, which we demonstrate analytically, and conductivity calculations show a hard gap with the 
current dual to the active gauge field superconducting, while the other is insulating. In the other 
limit, conductivities reveal a soft gap with the roles of the gauge fields reversed. 
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1 Introduction 



States of ordinary metals are well-described by Fermi liquid theory, consisting of weakly-coupled, 
long-lived quasiparticle excitations around a Fermi surface. It is of interest, however, to move 
beyond standard Fermi liquid theory; indeed, certain strongly correlated electron systems, such as 
the "strange metals" arising in high- Tc cuprates [1, 2J and heavy fermion systems [3j, are known to 
possess a Fermi surface from photoemission experiments, but the associated gapless excitations are 
not long-lived. Such systems, with Fermi surfaces but without ordinary quasiparticles, have been 
called "non- Fermi liquids" , and developing frameworks to characterize them is important from the 
theoretical point of view. 

The gauge-gravity correspondence [H EJ |6] is a promising avenue to investigate these phenom- 
ena, as strongly-coupled systems which are difficult to characterize in non-gravitational variables 
are described holographically by black hole geometries that may be far more tractable calculation- 
ally. Holographic Fermi surfaces were first investigated in [71|8l|9], and have been the subject of 
substantial study. Non-Fermi liquid systems have indeed been obtained; a study of Fermi surface 
behavior for general dimension, mass and charges was carried out in |10j . 

Much of the work so far has followed a "bottom-up" approach: a convenient effective Lagrangian 
in the gravity theory is postulated, and its consequences worked out. However, while such an 
approach is quite valuable, it does not make direct contract with supergravity or string theory. 
The potential drawbacks are twofold. First, the exact nature of the dual field theory remains 
unknown, as only the symmetries can be reliably established without the string theory map, so one 
is ignorant of exactly what system is supporting the Fermi surface. Second, without the guidance 
of a known string theory construction, one may worry that the masses or charges selected may 
contain a hidden instability or other secret unphysical property, and therefore that the results may 
be less trustworthy. Therefore, a "top-down" approach, where one starts from a known string 
or supergravity construction, is naturally valuable both for understanding the system and having 
confidence in its validity. Previous top-down approaches to Fermi surfaces include |1H [T2] , which 
studied fermions realized on probe branes, and |13 1 ll4 t [T5]. which studied the gravitino in the gravity 
multiplet of minimal supergravity and found no Fermi surface singularity; [16] (see also |17l I18j ) 
discusses rotating black holes with zero entropy in supergravity and relates them to Fermi surfaces 
from a different perspective. For other studies of finite density systems using the gauge-gravity 
correspondence, see for example [l9]-[l5]; recent reviews appear in [36 l WI \ H8]. 

Many known holographic examples also have a nonzero entropy remaining even at zero temper- 
ature. This characteristic is somewhat mysterious physically and is unlike most strongly correlated 
systems that one might want to apply the gauge/gravity correspondence to. Thus, it is also useful 
to pursue examples of holographic Fermi surfaces where the zero-temperature entropy vanishes. 

In [l9], we presented top-down constructions of Fermi surfaces in five- and four-dimensional 
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maximally supersymmetric gauged supergravity, dual to four- dimensional = 4 Super- Yang- Mills 
(SYM) theory and three-dimensional ABJM theory, respectively. The work [49j studied for each 
case a single mode in a single background corresponding to all chemical potentials equal, and was 
primarily concerned with establishing the existence of top-down holographic Fermi surfaces. It is 
natural, however, to wish to go beyond a single case to study a wide variety of fermionic modes, 
in a variety of backgrounds with varying chemical potential, and attempt to classify the diversity 
of Fermi surface behavior. Moreover, the explicit knowledge of the dual field theory in top-down 
constructions means an understanding of the dual zero-temperature states is within reach, and 
understanding a broader class of solutions can only help make this identification. In this paper, we 
undertake such a study for the case of A/" = 4 Super- Yang- Mills. 

= 4 Super- Yang-Mills has an SO (6) R-symmetry within which three distinct chemical po- 
tentials may be turned on. Here we study the class of extremal black brane geometries dual to 
zero-temperature states in which two of these chemical potentials are set equal to one another, 
while the third may vary separately; the dimensionless ratio provides us with a one-parameter 
family of geometries. The cases where one chemical potential or the other is set to zero are also 
of interest, and we study extremal limits of these as well: one is a non-thermodynamic BPS state 
on the Coulomb branch, and the other is a black hole with vanishing entropy at zero temperature. 
We study all the spin- 1/2 fermionic modes of five-dimensional maximal gauged supergravity — 
corresponding to the smallest-dimension fermionic operators in the large- A^, large-coupling limit — 
that do not mix with the spin-3/2 gravitino fields. There are 18 such modes, along with their con- 
jugates, obeying 8 distinct Dirac equations. Characteristic of these equations are scalar-dependent 
(and therefore position-dependent) mass and Pauli terms, absent in the usual class of bottom-up 
models. 

The extremal geometries with two distinct chemical potentials, called the 2+1-charge black 
holes, are of a class with a double pole at the horizon studied in general in In this paper 

we numerically solve the various Dirac equations in these backgrounds, looking for Fermi surface 
singularities. We find that fermions dual to Tr F^X operators made from a field strength and a 
gaugino have no Fermi surfaces, while fermions dual to Tr AX made from a gaugino and a scalar may 
have two, one or zero Fermi surfaces, depending on the net charge. Double fermi surfaces appear as 
nested spheres in /c-space. These systems all behave as non-Fermi liquids, with excitations that are 
not stable in the low-energy limit and hence are not true quasiparticles; we determine dispersion 
relations and width per energy ratios for these excitations. The fermions also generically possess 
so-called oscillatory regions manifesting log-periodic behavior where the Fermi surfaces disappear 
[HllTO], and we identify these as well. 

There are two interesting limits of our family associated to vanishing zero-temperature entropy. 
One is associated to taking the two identical chemical potentials much larger than the third, and 
we show that in this limit one of our Fermi surface singularities approaches a marginal Fermi liquid 
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(MFL), which occupies the border between Fermi and non- Fermi hquids, and was postulated to 
describe the optimally doped cuprates [1]. We are able to show analytically that this behavior is 
exact, and is controlled by the extremal limit of the "1-charge" black hole where the two equal 
chemical potentials strictly vanish, which is a non-thermodynamic BPS Coulomb branch state |50j . 
We solve the fermion fluctuation spectrum in this background, completing the analysis of [51]. We 
calculate conductivities in this geometry as well, and find a hard gap manifested by a step function, 
and argue that the current dual to the active gauge field is superconducting, while the other is 
insulating. We also determine the conductivities in the 2+1-charge backgrounds near the marginal 
Fermi liquid limit, which we may think of as "doping" the 1-charge background with a little of 
the other charge, and see that the delta function in the conductivity associated to superconducting 
behavior in the 1-charge background broadens into a Drude peak. 

The extremal limit of the "2-charge" black hole where the identical chemical potentials are 
nonzero but the other is zero is also of interest, for it is also dual to a system with zero entropy 
at zero temperature. We study the fermion fluctuation equations in this geometry as well, and 
find that while two modes have Fermi surface singularities, the form of the nearby small energy 
fluctuations is damped rather than infalling. This system displays a soft gap with a rapid but 
smooth rise in conductivity, again with a zero-frequency delta function in the active gauge field, 
although this may be a result of translation invariance rather than superconductivity. 

In section [2] we review the charged black hole solutions we will use, as well as their field the- 
ory duals, and in section |3] we show how these solutions are embedded within five-dimensional 
maximally supersymmetric gauged supergravity. In section |4] we determine the Dirac equations in 
these backgrounds for the spin- 1/2 fields of gauged supergravity, as well as describing their dual 
operators. In section [5] we discuss the near-horizon behavior of the Dirac equations and make gen- 
eral statements about oscillatory regions and Fermi surfaces, and in section [6] we solve these Dirac 
equations in the 2+1-charge backgrounds and obtain Fermi surfaces, whose properties we plot as 
functions of the ratio of chemical potentials. Section [7| studies the 1-charge black hole, enabling us 
to analytically derive the marginal Fermi liquid limit, as well as calculating the conductivities in the 
1-charge and nearby 2+1-charge backgrounds. In section[8]we determine Fermi surface singularities 
and conductivities in the 2-charge background, and in section |9] we conclude. 

2 Charged black hole solutions 

N = A Super- Yang-Mills theory with SU{N) gauge group is dual to type IIB string theory on 
AdS^ X S^. In the limit of large N and large 't Hooft coupling, string theory may be approximated 
by type IIB supergravity. The Kaluza-Klein reduction of type IIB supergravity on the five-sphere 
may be consistently truncated to only the lowest-mass modes constituting a single five-dimensional 
multiplet, resulting in five-dimensional maximally supersymmetric gauged supergravity. Each mode 



3 



in the gauged supergravity theory is dual to a known operator in the field theory. 

The reduction on the leads to an 50(6) gauge group in five-dimensions. Black brane so- 
lutions are known that are charged under the SO{2) x SO{2) x SO{2) = U{l)a x U{l)b x U{l)c 
Cartan subgroup. (We will use the terms "black brane" and "black hole" interchangably; all our 
geometries have three-dimensional spatial translation invariance.) The gauged supergravity theory 
also contains 42 scalar fields; general black holes with three unequal charges (Qa^Qb^Qc) source 
two distinct scalars and leave only the Cartan subgroup unbroken. The consistent truncation of 
the full theory to include the metric, the three gauge fields and the two scalars is called the STU 
model f52]. 

A simplification of this situation, which we pursue here, is to take two of the three charges 
equal; we take Qb = Qc and define Qi = Qa, Q2 = Qb = Qc- In this case we have a larger 
preserved symmetry, SO{6) — t- S0{2) x 50(4), and only one scalar field is turned on. The effective 
Lagrangian for the metric, two gauge fields and single scalar (in mostly minus signature) is 

e-i£ = -R+ 1(50)2 + + ±.^=^ _ e^f^.r'^ - 2e^^F^,F^'' - 2^^^^^'^^ f^.Fp^A^ . (1) 

The unconventional gauge field normalizations will be convenient for matching to maximal gauged 
supergravity. Here we will study a class of black brane solutions 

ds^ = e^^'^'\h{r)dt' - dx^) - ---dr^ , 

h{r) (2) 

a^dx^^ = ^i{r)dt, A^dx^ = $2 [r) dt , cj) = (p{r) . 

Particular limits that have been studied and will play central roles here include the "1-charge black 
hole" (Q2 = 0), the "2-charge black hole" (Qi = 0) and the "3-charge black hole" (Qi = Q2). 
We will refer to a solution with general Qi, Q2 as a 2-|- 1-charge black hole. We now review these 
solutions. 
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2.1 2+1-charge black holes and 3-charge case 

The 2+1-charge black hole involves both gauge fields and as well as the active scalar. It 
satisfies the ansatz ^ with solution 

A{r) = log^ + ilog fl + ^log fl + 



L 6 \ y 3 \ 

B{r) = -log^ - ilog fl + |log fl + 



L 3 \ y 3 \ 



where we have traded in a mass parameter for the horizon radius rn, which is the largest solution 
of h{rH) = 0. In the limit of large r, the solutions approach anti-de Sitter space with radius L. 
The temperature and chemical potentials are determined by the usual methods to baj 



-t = / > /^l = / , ^2 = ^2 ' 'y^) 



and the entropy and charge densities are 

An extremal black hole exists for any solution setting the temperature to zero with chemical po- 
tentials nonzero. This requires a non- negative solution to 

rl = -jQ\ + ^QlQl-\Ql, (6) 

and it is easy to see for nonzero Qi and Q2 such a solution always exists. Thus for general charges 
in the 2+1-charge black hole there is an extremal limit. The entropy density is generally nonzero 
at extremality. In practice it will be more convenient to solve the extremality condition ^ by 
eliminating (52, 

Ql = ''^ + rl, (7) 

so the horizon radius rn remains an explicit parameter. 

Perhaps counterintuitively, the ratio of chemical potentials /x/j = ^1/112 at extremality decreases 



^The chemical potentials are defined relative to canonically normalized gauge fields. 
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as Q1/Q2 (or Qi/rn) increases, 



M2 



while the ratio of charge densities 

Pi Qi 



(9) 



P2 V2Q2 ' 

more expectedly moves in the other direction. Furthermore, we see the expression ([s]) is bounded 
above at /x/j = 1; no extremal black holes in this class exist with larger ratio. The Q1/Q2 — ^ and 
Q1/Q2 ^ 00 (or Qi/rn — )• and Qi/rn — )• 00) limits of the extremal solutions are related to the 
1-charge and 2-charge black holes, as we describe momentarily. 

The 3-charge black hole is the equal-charge special case of a 2+ 1-charge solution: 

Qi = Q2 = Q. (10) 

which holds when 

A^ = a^ $2 = ^1 ■ (11) 

The scalar field thus vanishes for the 3-charge solution, which solves the simplfied effective action 



e-'C = R+^- ^UuF" - 2e^-P^^U,f,,ar . (12) 



The 3-charge black hole has an extremal limit as well, with rjj = Q /2, corresponding to 

Unlike the 2-charge and 1-charge black holes we describe next, the 3-charge geometry behaves as 
just another member of the general class of 2+ 1-charge black holes. 

2.2 2-charge black hole 

The 2-charge black hole comes from setting Qi = in the general 2-|-l-charge solution ([s]); this 
removes the gauge field from the system. The solution is 



(14) 



^(r)=log^ + ^log(l + ^) , i?(r) = -log^-^log(l + ^ 
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The thermodynamics is 



rH _ V2Q2 _ '^H , 2 , na\ _ ^Q^s 



^-vrL^' ^^ = ^' '=^.^^H + Ql), ^2 = ^^, (15) 

with Hi = pi = 0. There is an extremal solution for rn = 0; most interestingly for us, we notice 
that in this limit the entropy vanishes, 

ST=o = 0. (16) 

For the extremal solution the scalar diverges at the horizon, suggesting new light modes appear and 
the geometry is strictly IR incomplete; however the missing modes are expected to be associated 
with unbroken abelian gauge symmetry based on a D3-brane configuration in ten dimensions |531 l50j 
(where the scalar describes the squashing of the S^) which are 0{N) and so subleading. 

2.3 1-charge black hole 

The 1-charge black hole, conversely, comes from setting Q2 = 0, which removes the gauge field 
from the problem. The solution is 

^(r) = log- + -logh + ^j , i5(r) = -log---log(l + ^ 



h{r) 



(17) 



The temperature and remaining chemical potential become 



(18) 



with //2 = 0, and the entropy and nonzero charge density are 

Note there is no longer any solution for an extremal black hole. Instead, the limit r// — t- of the 
1-charge black hole is a BPS configuration with no gauge fields and no horizon function, preserving 
four-dimensional Lorentz invariance with only the scalar running. It is actually the so-called n = 2 
renormalization group flow solution studied in |50j . corresponding to a zero temperature and zero 
chemical potential state on the Coulomb branch of = 4 Super- Yang-Mills; our scalar is related 
to the one in |OTj by = — 2/Li, the parameter £ of that solution is simply Qi and the variable v 
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used there is given by 



(20) 



r2 + Qf • 

The 1-charge black hole is also a Poincare patch limit of superstars considered in [53]. Again the 
scalar diverges and the geometry is IR incomplete, but the lift is known to be a disc of D3-branes 
|50j and the missing modes should again be 0{N) in number. 



2.4 Relation of extremal geometries 

Up to overall rescaling of coordinates, the extremal 2+1-charge black holes are a one-parameter 
family of solutions determined by Qi/vh (or equivalently Q1/Q2 or Q2l'ru). It would seem natural 
that the extreme limits of the one-parameter family are the extremal 2-charge black hole on one 
side and the Coulomb branch solution on the other. In fact, order of limits issues make these 
identifications somewhat subtle. 

The Coulomb branch solution follows from taking Q2 ^ ^ strictly to reach the 1-charge black 
hole, without first requiring extremality, and only afterwards taking ru — t- 0. Thus the ratios of 
the parameters are 

§^^0, ^^0, — ^0 (Coulomb branch) . (21) 

Q\ Qi ru 

On the other hand, if we first take the extremal limit and afterward take Q2 ^ ^ with Qi held 
finite, then Q2 and th are linked by the extremality condition and we have 

§^^0, ^^0, — ^1 (extremal 2 + 1) . (22) 

Ql Ql TH 

The same ratios result from Qi — )■ 00 with th fixed when one has eliminated Q2 for extremality 
using Q. The Coulomb branch solution also has /X2 = 0, while the limit of the 2+1 black holes 
has instead ^i//U2 — )• 0. 

Thus the limits are not the same. As it happens, despite the difference in the limits, the 
metric, gauge field <I>i and scalar of the 2+1-charge extremal black holes nonetheless approach the 
Coulomb branch solution in the r// — )• limit. The only field that does not is $21 which instead of 
approaching <I>2 — ^ as one should have for the 1-charge black hole, instead achieves the nonzero 
value 

$2 ^ ^ . (23) 

Thus a general field in the rj^ — t- limit of the 2+1-charge extremal black holes will not explore 
the Coulomb branch geometry. A field that is unaware of <I>2, however, will find its dynamics 
approaching that of the Coulomb branch. We will see this come in handy later. 

There is a similar story for the 2-charge extremal black hole. The true 2-charge extremal case 
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involves Qi — )• strictly, followed hy rn ^ with Q2 finite. This implies 

^^0, "Jf^O, (extremal 2) . (24) 

Q2 Q2 ru 



Indeed we have in this case 



Q1Q2 



(extremal 2) , (25) 

'H 

for any a and /?. In addition the chemical potentials are /^i = 0, /i2 7^ 0. 

On the other hand, if we first take the extremal limit of the 2+1-charge black holes and then 



take Qi — )• with finite, we find Q2 — ?• 00. Although the relations (24) all follow in this case, 
we also have 

^ V2 (extremal 2 + 1) , (26) 

whereas in the true 2-charge black hole, QiQ2/f'\j = 0. In addition the chemical potentials approach 
/^i//^2 1- Thus again the limits are not the same. 

Examining the solutions, again only one field fails to agree in the two limits: this time it is <I>i, 



which is zero in the true 2-charge black hole, but in the limit described around (26) approaches 

(27) 

Any fermion insensitive to the gauge field would thus not notice the difference in limits; however, 
unlike the previous case, we will have no such fermions. 



2.5 Field theory duals 

The dual field theory to the extremal charged black hole solutions is = 4 Super- Yang-Mills with 
gauge group SU {N) at large number of colors and large 't Hooft coupling, at zero temperature and 
nonzero chemical potentials /^i and jjL2- 

The field content for AA = 4 SYM is a gauge field, adjoint Majorana fermions A^, a = 1 ... 4 in 
the 4 of 50(6), and adjoint scalars X in the 6 of 5*0(6). The charges of the dual gauginos and 
dual scalars Zj = X2j~i + iX2j, j = 1,2,3 under the Cartan subgroup U{l)a x U{l)b x C/(l)c are 





Ai 


A2 


A3 


A4 


Zi 


Z2 


Z3 


qa 


1 

2 


1 
2 


1 
2 


1 
2 


1 








Qb 


1 
2 


1 
2 


1 
2 


1 

2 





1 





Qc 


1 
2 


1 
2 


1 
2 


1 
2 








1 



(28) 



9 



The charges qi = qa and q2 = Qb + Qc associated to the chemical potentials and /i2 are then 





Ai 


As 


A3 


A4 


Zi 


^2 


Zs 




1 

2 


1 

2 


1 
2 


1 

2 


1 










1 


-1 











1 


1 


93 


3 
2 


1 
2 


1 

2 


1 

2 


1 


1 


1 



(29) 



where we also included qs = qi + q2- 

The nonzero scalar field (f>{r) in the gravity background Q leads to a vacuum expectation value 
for the dual operator in the 20', 



O 



20' 



Tr(-2|Zi|2 + |Z2p + |Z3| 



(30) 



implying expectation values for the individual scalars. The expectation value for O20' vanishes at 
the 3-charge black hole fiji = 1/^/2, and has opposite signs on either side of this point. 

We will introduce the fermionic fields we study in this paper and the operators they are dual 



to, including their charges, in section 4.2 



3 Charged black holes from maximal gauged supergravity 

In order to derive the fermion fluctuation equations in the black brane geometries, we first obtain 
the embedding for the Lagrangian ([T]) in the full maximal gauged supergravity theory. This section 
is somewhat technical and may be skipped by a reader anxious to get to the results. 

3.1 Scalar coset representative 

The scalar fields of Z) = 5, A/" = 8 gauged supergravity [55] parameterize an E^f^Q-^/U Sp{S) coset 
manifold, and the SO{&) gauge group is embedded in both the numerator and denominator groups. 
Under 5*0(6) the scalars transform in the representations 20' ® 10 © 10 © 1 © 1. The -E'6{6) is 
conveniently parameterized in terms of an SL{6, M) x SL{2, M) subgroup, with SO{6) embedded in 
(5*1/(6, M) in the natural way. Elements of the 20' sit in the coset SL{6,M)/ S0{6), while the singlet 
complex dilaton sits in the SL{2,M) / SO{2) and the 10 © 10 sit in the off-diagonal blocks. 

The scalar arising in the 2+1-charge black holes lives in the 20' and may be described in terms 
of the SL{6,M) element 

</> —'!> —0 -0 — /oi\ 

S = diag{e^,e^^,e2v^,e2v^,e2^/6,e2V6} ^ (31) 

while the SL{2, M) element is S" = I and the off-diagonal terms are zero; (pix) is the dynamical 
scalar field. Decomposing the S*L(6, M) indices /, J = 1 . . . 6 into x,y = 1,2 and i,j = 3, ... 6 and 
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denoting the SL(2, M) index a, /3, we then have 



sy = e^^5y, 5/ 



'2^6 J-? 



(50/ = 



(32) 



The -£^6(6) elements U are defined a^ 



u 



IJKc 



Uk. 



alj 



uJ^ = S/iS')J, (33) 



and the scalar coset representatives are then 



V 



IJab 



V, 



ab 



la 



4^/2 



abjT 

a 



(34) 



Here a, 6 = 1 ... 8 are USp{8) indices raised and lowered with the symplectic metric Qab and its 
inverse Q""^: 

Aa = ^abA" , A'' = , (35) 

where $7^^ = 17"^ = — fij^ The Fx, X = 0, 1, ... 6 are hermitian, antisymmetric and pure imaginary 
8x8 Euclidean F-matrices obeying 



{Fx,Fy} = 25xY , 



(36) 



and related to each other by 



Fo = iFi • • • Fe . 



(37) 



The F-matrices act on the USp{8) space as Clebsch-Gordan coefficients. Following [55], we separate 
out F/, / = 1 . . . 6 from Fq, and define 



Tij = ^{TiTj - FjFj) , F,„ = (F,, iF.Fo) , 



(38) 



where a = 1,2. It is straightforward to show that F/j and F/q are also antisymmetric as 8 x 8 
matrices. The USp{8) metric Qab may be written as 



ab 



(39) 



^By using the inverse of S in the first term we correct a typo in equation (5.34) of [55) . 

counterintuitive fact is that is not what you get when you raise both indices of Q,ab according to this rule; 

you get n'"'n'"^ncd = n''". 
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The coset representatives (34) then become 



^xa ~ 2^/2 ' ~ 2\/2 v-"- «oJ ) 



(40) 



and the inverses are 



(41) 



2x/2 ^ ' 2^/2 

3.2 Composite connection and scalar kinetic terms 

The covariant derivative of D = 5, A/" = 8 gauged supergravity acts on both USp{8) indices with a 
composite connection Qf^^J', and on 5*0(6) indices with the fundamental gauge fields A^jj = A^^jjy. 

D^Xal = d^Xal + QnaXbi - gA / Xaj , 

(42) 

D^X'^' = d^X-' - Q^.'^X"' - gA/jX-' . 



The covariant derivative of the coset representative defines the quantity P 



abed 



p^abcd ^ yabABjj^y^cd ^ (43) 

where AB is a formal index pair in the 27 of -56(6) ' whose contraction is defined as 

V-'^^'D^VABab = V'^'jjD^V'i, + V-'^'^D^Vi^ab , (44) 



plectic traceless, denoted p "■'"^ = p^"' this leads to 



and where the composite connection is defined so that P^"-^'^'^ is totally antisymmetric and sym- 

[abcd\ I _ 



In our background the composite connection is 

b _ 9 „ ,^ .ab 9 



(45) 



Q,a " = -T^a.iFur - |^/.(r34 + Tse)"" ■ (46) 



2 2 

The scalar kinetic terms in the gauged supergravity theory are 



ll...X.r=Y^P^abcdP^'''''\ (47) 
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(48) 



and p ^^^'^ evaluates to 
leading to the kinetic term for the neutral scalar (p, 

/:scalar = J(5</.)2. (49) 
o 

As we will see, this unconventional normalization will become canonical when we rescale the gauged 
supergravity Lagrangian to match ([T]). 

3.3 Scalar potential 

The gravity /scalar sector is then given by the Lagrangian 

The potential may be determined in terms of a few U Sp{8) tensors that depend on the scalar; these 
tensors will also enter into the fermion couplings. These start with the tensor T'^^^^, 

= (2y^^- VheJK - K7„"^ V,, VaiJ , (51) 
which is evaluated to be (with one index lowered), 

Tabcd = ^evl (ror,,)'^'(r/,r' + ^evs (ror,,)'^''(r,,)^'^ . (52) 

lb lb 

It is easy to see that this is symplectic traceless in the first pair of indices and in the second pair, 

^""^Tabcd = ^'^'^Tabcd = , (53) 

as well as antisymmetric in the second pair but symmetric in the first, 

T(ab)cd = Tab[cd] = Tabcd ■ (54) 

The tensor Tab is defined as 

Tab^T'^abc, (55) 
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and evaluates to be 

/ 15 -'^ 15 20 \ 
Tab=[-^eV-e--eV-.j5ab, (56) 

and finally the tensor Ajjcd is defined as the antisymmetrization and trace-removal over the final 
three indices of Tabcd, 

Aabcd = Talpcd]\ , (57) 

which can be shown with a little calculation to be 

^abcd = 2 {Tabcd + Tacdb + Tadbc) + g {Tabbed + Tac^db + Tad^bc) ■ (58) 

It is also straightforward to see that T"^^'^, T"^ and A'''""^ have the same values for their components 
and Aabcd- The five-dimensional M = 8 gauged supergravity potential is then defined 
in all generality as 

V=^ (^AabcdA'^"''' - ^TabT'^') , (59) 

and is in our case 



F(</,) = -^eVi - . (60) 



2 4 

The superpotential VF((;^) may be defined as Wab = W{(p)6ab, where 

Wab = ^Tab , (61) 

15 

becoming in our case 

-0 1 M 

il/(^) = -eV6 - -ev^ . (62) 
The potential may be recovered from the superpotential by the usual gauged SUGRA formula 

3.4 Gauge fields 

Turning to the gauge fields, the gauge kinetic term coming from J\f = 8 SUGRA is 

6 -^gaugekin — '^^ixi/ab^^ j (64) 

where we have set the two-form fields to zero, and 

F^,''" = F^.^uV'^'K (65) 
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Using our expressions (40) for the coset representatives, we obtain for the three Cartan gauge fields 
Af_ii2, A^34, Afj,5e, 

1 -40 ^ 'i'i> 1 20 

e-^/:gaugekin = - F^^uF^'^^^ - -e v5F^,34F^-34 _ ^eVSF^,56F^-56 _ (66) 

Let us define 

afi = Afj^u , An = A^u = ^^56 , (67) 

and then we have 

e-^/:gaugckin = -\e^fn,f'''' - hv^F^.F'^'' . (68) 
There is also a Chern-Simons term, 

Ccs = -^e^'""'^e'-''''^'''''Fijn,FKLpaAMNr , (69) 
where we have dropped 0{FA^) and 0{A^) terms which vanish for us. This evaluates to 



Ccs = -^-e^'P'^- {2UuFp.Ar + F^^Fp^ar) = -^e'^'""'^ f^.F^^A^ , (70) 



where in the second equality we integrated by parts and dropped a surface term. The complete 
gravity /scalar /gauge Lagrangian is thus 

1 1 1 n o / I 1 -2'^'\ 1 -40 1 20 1 

(71) 

If we identify the gauged supergravity coupling g with the AdS radius L, 

(72) 

the Lagrangian precisely matches ([T]) up to an overall factor 1/4, which can be absorbed into a 
definition of the gravitational constant. 



4 Fermionic fluctuation equations 

We now study the Lagrangian for spin-1/2 fields in = 8 gauged supergravity, and extract 
the Dirac equations for fermionic fluctuations in the backgrounds of the black holes discussed in 
section [2] A reader wishing to skip the details of gauged supergravity may proceed to the general 



Dirac equation (80). 
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4.1 Spin-1/2 Dirac equations 

The spin-1/2 fields are the 48 x"'^^) which are both antisymmetric and symplectic traceless, = 
0. The 48 fermions are complex, but not independent; they are related by the symplectic Majorana 
condition: 

^abc^C^-abc-^T _ (73) 

where C is the conjugation matrix defined via 

(7^)^ = C-i^'C-^. The barred version of the fermi 

field is defined via conjugation as 

r'^=(Xa6c)V- (74) 

Note this definition relates a conjugate fermion with indices up to the original fermion with indices 
down. Thus x"^^ and Xafeo which are a priori independent fields thanks to raising and lowering by 
the symplectic matrix (for example, x^^^ = X123 7^ X567) are related as one another's conjugates 
by the symplectic Majorana condition. Hence for the 48 spin-1/2 fields, only 24 are independent; 
the remaining are related to the conjugates of the first 24. Completely analogous relations hold for 
the 8 gravitini ?/^^. 

The kinetic term for the fermi fields is 

e-^C^n=^r'"l^D^Xabc: (75) 

where the covariant derivative acting on the x"'^'^ is 

D^Xabc = ^/iXabc + Q^a Xdbc + Q ^b Xadc + Q i_ic Xabd , (76) 



where contains the spin connection and Q is the composite connection ( 46 ) ; this provides 



the gauge couplings for the fermions. The spin-1/2 fields also have mass terms of the form, 

e~^/^mass = -^X"''"^ (^2^bcde " -^^bdTce^ Xo!^^ ) (77) 



where Tab (56) and Aabcd (58) are tensors on the scalar space defined previously, as well as Pauli 
terms 

e-^£pauli = iF^.'^'Xacdl^^Xb"" 

(78) 
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e^^.(ri2)"' + e^F^,(r34 + Tse)'^') XacAl^^'Xb'' ■ 



We are interested in spin-1/2 fields that do not couple to the gravitini; we will discuss how to assure 
this momentarily. For such fields the quadratic Lagrangian in a 2-|-l-charge black hole background 
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takes the form 



2 



X: 



(79) 



where the vector on x stands for the 48 fields encoded in the abc indices. One may diagonahze and 
normalize the kinetic term K to the identity. The eigenvalues of M, Qj and P, may then be found. 
The equation for a fermion eigenvector is then 

( -ti> 24, \ -24, 4> 

mie^ +77126^]+ gqij^Qf, + 9927''^^ + ipie ^ Uul^"" + ip2e ^ F,,^-/'"' x = . 

(80) 

Here we have extracted a factor of g from the eigenvalues of Qj and exponentials depending on the 
scalar from the eigenvalues of Pj, as well as writing the scalar-dependence of the M eigenvalues, 
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m = g \ rriie ^ + m2e ^ 



The fermion is then characterized by the six rational numbers {rrii, qi,Pi 



(81) 



It is straightforward to show that if a field Xabc obeys (80) with couplings {rrii, qi,pi), then 
the conjugate field/Majorana partner C{x"'^'^)'^ = x"^^ obeys (80) with couplings {—nii^—qi^pi)] 
the eigenvectors thus come in pairs. A negative value of m corresponds to a fermi field whose 
field theory dual has opposite chirality. In practice, one may work only with positive masses, by 
changing the Clifford basis of the negative mass modes, flipping the sign of all 7-matrices. This 
effectively changes the sign of m and pi. Thus instead of {—rrii, —qi,Pi) one may solve the equation 
for the conjugate/partner field with {rrii, —qi, —Pi), bearing in mind that the chirality of the dual 
field is opposite. 



4.2 Dual operators and "maximal" spin-1/2 eigenvectors 

The 48 spin-1/2 fields transform in the 20 -|- 20 -|- 4 -|- 4 representations of S0{6). The fields in the 
20 are dual to A/" = 4 Super- Yang-Mills operators of the form, 

X"^" ~ Tr (AX) . (82) 

where A is the gaugino in the 4 and X is the adjoint scalar in the 6. There are twelve such fields 
with the mass values (mi, 1712) = (1/2,-1/4), four with mass values (—1/2,3/4) and four with 
mass values (—1/6, 5/12); for all these cases the mass with vanishing scalar becomes 

m2o(0 = O) = | = ^. (83) 
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Meanwhile the spin- 1/2 fields in the 4 are dual to the operators 



X ~ TV (F+A) , (84) 

where is the self-dual part of the SU (N) field strength. These fields have mass values (m-i, 7712) = 
(1/2, 1/4); this mass as a function of the scalar fields is proportional to the superpotential and for 
vanishing scalar becomes 

m4(0 = O) = | = A. (85) 

We would like to study spin- 1/2 fields that do not mix in the quadratic Lagrangian with the 8 spin- 
3/2 gravitino fields ip^, which transform in the 4 + 4. A number of such fields can be determined 
by exploiting the conservation of the preserved U(l) charges. Since the metric and the scalar (p are 
neutral, the quadratic fermion action can only couple fermi fields with equal and opposite charges 
under any preserved U(l)s. So if we identify a field x"'^'^ that has a charge possessed by no we 
can be assured that x"'^'^ does not mix with the gravitino in the quadratic action. 

We can find the most decoupled spin-1/2 fields by considering a more general background: the 
U{l)a X U {l)b X U{l)c solutions with A^i2, ^^34 and ^^55 all nonzero and distinct. This background 
has two (neutral) scalar fields in general, but we will not need to consider the details of the solutions. 
In these backgrounds, each spin-1/2 field has three U{1) charges qa,qb, Qc and the above argument 
holds for all three. (In our background, qi = Qa and q2 = Qb + QcQ) 

The 4 decomposes into U (l)a x U (1)5 x U (l)c charge vectors (5,5,^), (2,-5,-5), (—5, 5, —5), 
and (—5, —5, 5), while the 4 is the negative of these; thus all three U(l) charges for the gravitini 
have magnitude 1/2. The 20, meanwhile, contains two copies of the weight vectors of the 4, and 
twelve distinct weight vectors where one of the three charges has magnitude 3/2. Thus of the 48 
spin-1/2 fields in 24 symplectic Majorana pairs, there are 12 pairs we can be assured do not mix 
with the gravitini due to having one "maximal" U{1) charge of magnitude 3/2. 



The charges of the dual gauginos Xa and dual scalars Zj are given in the table in section 2.5 A 
moment's examination reveals that the way to get an operator with one charge of magnitude 3/2 
is to make sure the nonzero charge of the scalar has the same sign as the corresponding charge in 
the gaugino. 

In the following table we list the 12 "maximal" modes x'^"'^^'^" with positive 53. All the maximal- 
charge fields have one of two sets of mass values. Some modes are right-handed; this is manifested 
when the eigenvalue of the mass operator comes out negative. For such fields we use the opposite- 
sign 7 basis as previously discussed, hipping the sign of the mass and Pauli couplings; these have 
already been flipped in the table. Such fields have right-handed duals and are indicated by x 
instead of x ™- the first column. 



*Note this is different from the rule for the black hole's charge, which is Q2 ~ Qb = Qc- 
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The fields ;!^(2'2'2)j t^(2'2'2) and ;^(2'2'2) have maximal 53 = 5/2, and thus in the background of 
the 3-charge black hole are the modes studied in [39]. 

The conjugates/symplectic Majorana partners of the modes in the table have the opposite 
chirality, and the signs of the Qi and pi flipped. 



4.3 "Overlapping" spin-1/2 eigenvectors 

The 24 remaining spin-1/2 fields have charges {qaQbQc) = (i^i =^5) i^)) each of these 8 charge 
assignments occurring with multiplicity 3; these modes overlap with the gravitini, which have the 
same 8 charge vectors. Although group theory does not rule out any of these fields from mixing 
with the gravitini, we can see explicitly that some do not. 

There are three types of x'^^'^ /^^i couplings in the gauged supergravity Lagrangian, a direct 
coupling mediated by the Aabcd tensor, a coupling with a scalar derivative 5^0 coming from P^abcd, 
and a Pauli-type term, which in our 2-|-l-charge background splits into two distinct couplings, one 
for each gauge field: 

We may diagonalize the mass matrix for each three-dimensional charge subsector of spin-1/2 fields. 
We find three distinct masses in each sector: one with mass values {mi, 1712) = (1/2, 1/4), one with 
(1/2, —1/4) and one with (—1/6, 5/12). Plugging the mass eigenvectors into the gravitino couplings 
( |87[ ), we find the (1/2,1/4) eigenvectors never couple to the gravitini; these are the higher-mass 
modes dual to Tr F+A. 

The (1/2, —1/4) eigenvectors in this sector have zero for most gravitino couplings, but two out 
of four of them notice the Pauli-type coupling. The (—1/6,5/12) eigenvectors, on the other 
hand, couple to the gravitino in almost every possible way; the only gravitino couplings they don't 
have are the two Pauli couplings just mentioned. 

In the table, the last four columns indicate whether the fermion in question has a gravitino 
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coupling of the sort indicated. Thus there are six distinct decoupled Dirac equations from this 
sector; the ones with mass values (1/2, —1/4), however, have the same Dirac equation as some of 
the "maximal" modes. Thus the only new Dirac equations we get from this sector are the four with 
mass (1/2, 1/4) which are dual to Tr (F+A), three of which are distinct. 

Again, the conjugate modes exist with opposite chirality and opposite signs for the qi and pi. 
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5 Solving Dirac equations, Fermi surfaces and oscillatory regions 

In order to look for holographic Fermi surfaces in A/" = 4 Super- Yang-Mills theory, we must analyze 



the various Dirac equations we have obtained. We do this in section 5.1 for a slightly more general 
background than we later require. Then we specialize to the 2 -|- 1-charge extremal black hole 
backgrounds, characterized by Qi/rn-, or equivalently, by = ^i//U2> and present the near- 
horizon asymptotic analysis followed by a numerical study that allows us to trace out the main 
features of the holographic Fermi surfaces as a function of /i^. In later sections we will consider 
more carefully the 1-charge and 2-charge limits. 



5.1 Solving the Dirac equation 

We now analyze the Dirac equation 



i7^V/, - m((/.) + gqi-i^'a^ + gq2l^Af, + ipie ^ f/.^j'"' + ip2e^F^^j'"' x = , (89) 
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in our class of backgrounds; we follow the basic conventions of |10j . though in opposite signature. 
The covariant derivative contains the spin connection, 



and we define 



X = e 



(90) 



(91) 



where for convenience we have placed the momentum k in the x-direction, and the e '^^h ^/'^ factor 
precisely cancels the effects of the spin connection. Choosing a 7-matrix basis including 



ifT3 



nr.)-'' 



i(T2 
, -ifT2, 



(92) 



we can define the projectors 



1 



1 



1 



-l)"z7^7*7^ 

with a = 1, 2, and characterize the four components of ^ as 

The Dirac equation decomposes into pairs relating to ^a- for each a: 
d,+ '^]m^_ = ^ [n(r) + {-lYk - v{r)] ^>^+ , 
5. - ^ ) = ^ [-n(r) + i-irk - v{r)] , 



(93) 



(94) 



(95) 



where we have definec2] 



v{r) = 2e~^{pie ^ dr*^! + p2e^ dr^2) ■ (96) 



We can decouple these into second-order equations, obtaining 



K± - F±K± + 



where 



me 
Vh 



F± = dr log 



^B-A 



(v{r) - (-l)°A:±u(r)) 



(97) 
(98) 



^One should not confuse v{r) with the variable v from [50! defined in (201 and used in section 7.2 
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One can see that (97) is invariant under 

Qi^ -Qi, Pi^ -Pi, uj ^ -uj , k ^ -k, (99) 

and so the equation of the conjugate fermion has the same solutions with {k,u}) — )■ {—k, —u). 

In the near-boundary hmit r — )• oo, where all our geometries approach anti-de Sitter space, the 
Dirac equation is dominated by the mass term. For \mL\ ^ 1/2, we have the solutions 

*«+~^a(A:)r'"^ + i?a(A:)r-"^-S ~ C„(A;)r'"^-i + Z)„(A:)r-'"^ . (100) 

The first order equations can be used to algebraically relate B to D, and A to C [56], 

_ LHu;+{-irk) LHu;-i-irk) 

- 2mL - 1 ' - 2mL + 1 ' 

For m > 0, A is the source term, and D the response; the case m < exchanges their roles. Thus 
as we have mentioned previously, negative mass corresponds to opposite chirality in the dual. 
For mL = 1/2, which is our primary case, we have instead 

^a- ~ Ca(fc)r-i/2 log ^ ^ D^ik)r-^/^ , (102) 
with the relation between B and D unchanged but the relation between A and C now 

Ca = L'^{oj + {-lTk)Aa. (103) 

The retarded Green's function G/j for the dual fermionic operator is defined in terms of the source 
A and the response D as 

Da = {GR)apA^ , (104) 
for a fiuctuation with infalling boundary conditions imposed at the horizon. Due to the decoupling 



of the Q components in (97), the Green's function is diagonal for us, with G22(w, k) = Gii{oj, —k). 



5.2 Near-horizon analysis: 2+1-charge case 

The extremal 2+1-charge black holes fit into a class of geometries with a double pole in the metric 
at the horizon, which has been studied in [lOJ. As r — )• r^/ we have the leading terms 

9ii^kQ, gtt^ -Toir -rnf , Orr ^ {L2f {r - r h)~'^ , 

(105) 
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where k^, tq, L2, Pi and (po are constants (we have eUminated Q2 to impose extremaUty using ([T])): 
ko 



'^CfT^^. ....... 



To = 


L 











2 ' 



2^Xffi' ^^"VS-^Q?;' (106) 



2rl ^ Qi^/2ri, + Qi 

Pi = / 7 P2 = — 



LQi^rjj + Qi 2LrH^rjj + Qi 

Let us consider the Dirac equation in the near-horizon Hmit at u = 0. We have 

gko , ^ ^ , 2kQ ( -^"^o ^ ^ ^ 

li {qxiix + g2P2) , V — — Pie ^ ^1 +p2e^^/52 

To T0-L2 V 

m —7- 'm{(j)Q) = g I mie ^ + ^ 



(107) 



leading to 

F± — . (108) 

r — rn 

The Dirac equation then simplifies to 

r — ru [r — rn) 

where i^fc is 



1 

dl^^^ + ^',± - ^^-ai , (109) 



with 



^k = \\[ + ( k/koY) {L2Y - 52 (giei + 9262)' (110) 



k = k- {-lY-^^ [ Pie ei +P2e^e2 ) . (Ill) 



^2 

and where we have defined the & 



I) 

ei = . (112) 

To 



taking the values 



' Qi(4r|, + Qf)' 4(44 + Qf) ' 



The solutions to ( 109 ) are then simply 

m ^ {r -th)^"^ . (114) 
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As described in |10j . these solutions possess a hidden near- horizon region approaching AdS2 x R^; 
appropriate boundary conditions are fixed by matching modes through this region. Defining the 
variables 

AL2 1 1 , , 

r-rH = 7, t=-T, (115) 

and taking the A — )■ limit with C, and r held fixed, we arrive at the metric 

ds^ = {-dr^ + d^) + kldx^ , (116) 

which is AdS2 x with AdS2 length L2, while the gauge fields become 

dx^" = ^dr , dx^" = ^dr , (117) 

identifying the as the coefficients of the gauge fields in the AdS2 limit. Imposing infalling 
boundary conditions in the AdS2 geometry and matching out into the outer region requires picking 



the ~ (r — th)'^'^'' boundary solution given above in (114). 



5.3 The oscillatory region, Fermi surfaces and excitations 



Consider the structure of the exponent z^^ (110). The effect of the Pauli terms is to shift the 
origin of the 3-momentum k to the quantity k defined in (111). The quantity is then a sum 
of positive-definite contributions from the near-horizon mass and the shifted momentum A;, and a 
negative-definite contribution from the effective near-horizon electric field coupling, 

(ge)efr = gqiei + 5^262 • (118) 

When the electric coupling is sufficiently strong certain values of k will produce an imaginary z^^. 
This is called the oscillatory region, associated with an instability toward the pair production of 
charged particles in the AdS2 region [57], and manifesting log-periodic behavior for the spinor 
excitations [HI fTO] . The boundary of the oscillatory region occurs at fcosc satisfying 

i^fcoec=0. (119) 

From (110), it is easy to see that 

~kL=['^)\{ie)l,-m'Ll) , (120) 

and thus the oscillatory region may only exist when the effective electric coupling {qe)'^g has greater 
magnitude than the effective mass rn^L^. It is thus significant that for the gauged supergravity 
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modes we study, the mass depends on the scalar field and hence is a function of /i/j; bottom- up 
constructions of fermion fluctuations witli a constant mass miss this aspect of the spectrum. (The 
effects of Pauli terms were considered in |58j-) As functions of /x/j, these take the form 



2 _ (V2gi/i| + g2(l -/^r))^ n 2 _ {■mijl - + m2^^l)'^ 

4(1-^|)(1+^|)2 ' "^^2- (,_^4^) 



while the proportionality factor is 

■fcoV_4r|,(2-/x|)(l + /x|) 



r , ,4 2n 2^ -2^1(1 + /^!), (122) 

where in the second equality we removed rfj- / in favor of the chemical potentials fii, fi2- Mean- 



while the shift between k and k ( 111 ) is opposite for the two components a = 1, 2 and is determined 
by the pi, 

L^HR \ ^-I^R (123) 

= (-l)"(2pi/ii+^/2p2;U2). 

In plotting our results for /cqsc and kp, it is useful to normalize these dimensionful parameters by 
a quantity meaningful in the field theory. The chemical potentials are a natural choice, and since 
< fii < fi2 for the extremal 2+1-charge black holes, the more convenient choice will be to plot 
the dimensionless ratio /c/^2 • The results then have the form 



k, 



osc 

^J'2 



A 



^osc I ^shift 
M2 / fJ'2 



(124) 



It is interesting to consider the expressions for k"^^^/ as ^r — >• 0: 




= ~ ""^ + (2m? - 4mim2 - ql) fx^ + V2qiq2 /^/j + • • • > (^^S) 



+fI^ + ^-4m,m2-^U.... (126) 



and as fiR — )• 1, 



\f^2 J 4(1 -^iii) V 2 V2 

We see that to leading order, the existence of the oscillatory region as ^r — ?■ is determined 
by a competition between q2 and mi, while the existence of the oscillatory region as ^ir — )• 1 is 
determined by a competition between qi and m2. In practice things are not quite so simple: we 
have g| = 4mf for more than half our fermions, and qf = Am^ holds for each and every one of them. 
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Thus the existence of the osciUatory region is generally determined by the subleading terms, which 
we therefore include above. The universal vanishing of qf — Am\ also means that the expressions 
for kosc are finite throughout the parameter space; this finiteness is another aspect that holds for 
our top-down fermions that would not be present in a generic bottom-up construction. 

Fermi momenta are defined as values of k = kp such that the source term in the near-boundary 



expansion (100) of the fermion vanishes: 

A{kF) = 0. (127) 

Fermi surfaces exist outside the oscillatory region. Varying /i/j can cause a Fermi momentum to 
pass inside the oscillatory region, at which point the Fermi surface ceases to exist. 

We will see that when a Fermi surface exists, the location of the Fermi momentum kp is often 
(but not always) right outside the oscillatory region. Thus while we do not have an analytic formula 



for the value of kp, we can gain some insight by using A^osc satisfying (124) as a proxy for the Fermi 
surface location. We will see there are several cases where kp ^ kosc holds quite well, but also a 
few where it does not. 

Once a Fermi surface is located, the properties of nearby excitations may be studied. Several 
of these depend solely on v^p ilQj. The retarded Green's function near the Fermi surface for small 
u! takes the form ^ 

GR{k, u) ~ , ^ i^,^ ,^ .2v^ , (128) 

with hi, h2 positive constants, k± = k — kp and the angle 7^ defined by 

7fc = arg (r{-2uk) (e-^--^ - e-2-('?-)off ) ) . (129) 

For non-Fermi liquids, which have u^p < 1/2, we can ignore the Fermi velocity vp term as sublead- 
ing. The dispersion relation between the excitation energy w^, and the momentum k± is then 

~ {k^y , (130) 

where the exponent is 

Furthermore the residue Z vanishes at the Fermi surface like 

Z^ik±y-\ (132) 
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Figure 1: Multiple fermi surface singularities correspond to nested spheres (left). When the signs 
of kp are opposite and the excitation spectrum is the same at both surfaces, one interpretation is 
a thick shell (right). 



Finally, the ratio of the excitation width to its energy is given by 



tan ( "^^^ 



tan ( TTZ 



A;_L > 0, 

fei < 0. 



(133) 



Different formulas hold for Fermi liquids, which have h'kp > 1/2; as we shall see, all our Fermi 
surfaces are non-Fermi liquids, with one interesting special case approaching the marginal Fermi 
liquid at u^p —7-1/2. 



6 Fermi surfaces in 2+1-charge black holes 

To obtain the locations kp of Fermi surfaces as a function of fiR = ^1/^2 for the various fermions, 
we numerically solve the decoupled second-order Dirac equation (97) at a; = 0, beginning at the 



horizon where we impose the positive sign exponent in ( 114 ) as a boundary condition, and searching 



for values of kp that cause the source term to vanish as (127). We arbitrarily solve for the spinor 
component a = 1; the other component a = 2 has identical results with k ^ —k. We plot the 
values of kp/fi2 vs. fiji for the Dirac equations with Fermi surfaces in the following figures, as well 
as i^kpj z and T/co for each case. 

A few general points before we consider each fermion in turn: 

• The fermions with asymptotic mass — ^ which sit in the 4 of 5*0(6) and are dual to the 
operators Tr -F+A, possess no Fermi surfaces. 
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For the fermions with asymptotic mass m ^ 7^, sitting in the 20 of S0{6) and dual to the 
operators Tr AX, the number of Fermi surfaces for each background can be classified by the 
total charge Qs = qi + q2- The two distinct equations with = 5/2 have two Fermi surfaces 
for at least part of the range of /i/j; the three distinct equations with q^ = 3/2 have one Fermi 
surface for at least part of the range of /i/j; and those with q-^ = 1/2 have none. This is in 
harmony with a general expectation that fermions with greater charge are more likely to form 
Fermi surfaces. 

The Fermi surfaces have < I'kp < 1/2, indicating they are all non- Fermi liquids. 

The Fermi surfaces for the most part exist close to the boundary of the oscillatory region; in 
each plot we include the values of kosc and shade the oscillatory region. As fiji is varied, kp 
may pass below A^osci at which point the Fermi surface ceases to exist. As the Fermi surface 
approaches the oscillatory region we have — t- 0, implying the dispersion relation exponent 
2; — )■ 00 as the system moves further into the deeply non- Fermi regime. 

At the 3-charge black hole point, denoted as 3QBH, the two Dirac equations q^ = 5/2 coincide, 
as do the three equations with 53 = 3/2. Thus the values of kp and /cqsc must match for 
fermions of the same type, as indeed they do. The qs = 5/2 modes at the 3-charge black hole 
point are the fluctuations studied in [49J, and we match the results found there; different radial 
coordinates and conventions for the sign of the gauge fields mean that Adhere — V^^there* 

When the spectrum of a fermion with charges {qi,Pi) contains a Fermi surface singularity at 
kp, the conjugate (antiparticle) equation with charges {—qi, —pi) will have a Fermi surface 
singularity at —kp (see equation ([99])); allowing k to point in the y and z directions will 
complete these antipodal points into a full three-dimensional Fermi sphere. For the examples 
where multiple values of kp are observed, this corresponds to nested Fermi spheres (see 
figure [TJ . 



• As noted in equation (133), one can examine excitations around Fermi surfaces with k± = 
k — kp > 0, corresponding to particles, and k± < 0, corresponding to holes; such excitations 
should be on the physical sheet of the complex w-plane jlOj. We found that excitations with 
A;_L < exist on the physical sheet for small u, but those with > do not, a manifestation 
of particle/hole asymmetry. In the following we plot T/co for k± < 0. The role of particles 
and holes is exchanged for the conjugate equation with charges {—qi, —pi)- 

We turn now to a discussion of the results for each of the five distinct Dirac equations for which 
we find Fermi surface singularities. 

Case 1: and x^^'^'^\ dual to Tr A1Z2 and Tr X1Z3 
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Figure 2: The values oikpj [x^ for case 1, Tr \\Z2 and Tr A1Z3. 




Figure 3: The values of f^^, z and F/o; for case 1, Tr A1Z2 and Tr XiZ^. 



These modes have (51,92) = (5,2), (^1,^2) = (i,0) and (mi, 7712) = (5,-j); kp is given in figure[2] 
and , z and T/uj in figure [3j 

The oscillatory region has the shape of a wavy band, as kosc/fJ'2 approaches a constant value in 
either limit, while /cghift approaches zero at ^ur — t- but a finite value at fiR — t- 1. This fermion is 
one of two with the largest total charge = 5/2, and has two Fermi surfaces for larger values of 
fiji; these track the oscillatory region boundary very closely. At the 3-charge point these fermions 
reduce to the results of [49J. The Fermi surfaces disappear into the oscillatory region at around 

^1R « 0.47. 

The values of are small throughout the range, resulting in a scaling exponent z for the 
excitations that never gets smaller than z > 5, and which (as it must) grows without bound as 
the Fermi momenta approach the oscillatory region. The ratio F/o; of excitation width also stays 
small, being bounded above by T/co ^ 1/10 and going to zero as /U/j — )• 1; thus in this limit the 
would-be quasiparticle excitations become more and more stable. 

In the full three-dimensional /c-space, these two Fermi surface singularities are completed into 
nested spheres. The excitation spectrum near both surfaces is the same; since the two kp solutions 
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Figure 4: The values of kF/fi2 for case 2, Tr XiZi. 




are of opposite sign, this imphes that the excitations inside the outer sphere match those outside 
the inner sphere. An interpretation of the two Fermi spheres is then that there is a "thick sheU" of 
occupied states in between them, with unoccupied states both on the outside and the inside of the 
pair (see Figure [T]). Put another way, there is a Fermi surface of antiparticles inside the surface of 
particles, canceling out in the overlap. It is also possible that the Fermi surfaces are associated to 
distinct would-be quasiparticles, but they must have opposite sign. 

Case 2: dual to Tr XiZi 

This mode has (91,92) = (§,1), {pi,P2) = ("i,^) and (mi, 777-2) = (-5,!); kp is given in figure [i] 
and ) -2 and F/cj in figure [sj 

As in case 1, both kosc and fcshift approach finite values at fiR — )• 1. At fiR — )• 0, however, the 
vanishing of 92 ~ ^777^ goes to zero, while kghiit reaches a finite value. Thus both 

branches approach kosc/fJ'2 — ^ l/\/2; the separation goes as 0{^'j^). 

This is the other fermion with largest total charge 93 = 5/2, and it also has two Fermi surfaces 
for much of the range of /i/j, again tracking the oscillatory region very closely. Unlike case 1, here 
the Fermi surfaces exist as far as can be determined for arbitrarily small fiR, while they disappear 
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into the oscillatory region on the right side of the plot, around ~ 0.93. At the 3-charge point we 
coincide with case 1 and the results of [l9]. The values of are again small, though the maximum 
is slightly greater than the previous case, and the exponent reaching a minimum around z > 3.5, 
with T/oj < 1/15. 

Interestingly, in this instance one of the Fermi momenta passes through zero, while the other 
does not. To the right of the zero, this case resembles case 1, where the two values of kp have 
opposite sign and one interpretation is an occupied shell in momentum space. To the left of the 
zero, however, the two Fermi singularities have the same sign of kp, implying the same excitations 
outside each Fermi surface; the nested spheres now have the same sign excitations. 

Case 3: and x^^hh^\ dual to Tr X2Z2 and Tr X2Z3 

These modes have (^1,^2) = (-5)2), (^1,^2) = (-jjO) and {1711,1712) = (5,-3); kp is given in 
figure [6] and u^p, z and F/w in figure [Tj 

Here the oscillatory region is a band as fiR — )• 0, with vanishing A:shift there due to p2 = 0, similar 
to case 1. Before can reach 1, however, the effective mass term dominates over the effective 
coupling and /cqsc ceases to exist. Hence for black holes with 0.85 < < 1, no oscillatory region 
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Figure 8: The values of /ci?//i2 for case 4, Tr A3Z3 and Tr A4Z2. 




Figure 9: The values of Vkp, z and F/a; for case 4, Tr A3Z3 and Tr \4,Z2. 



exists for any k. 

This fermion has total charge ^3 = 3/2, and there is one Fermi surface for most values of ^r. 
The Fermi surface emerges from the oscillatory region close to ~ 0.4, and above this kp tracks 
the value of fcosc for a while; however they soon separate and the fermi surface continues to exist 
even as — )• 1, once the oscillatory region has disappeared. Notice that kp goes through zero: 
this indicates that the Fermi surface shrinks to zero size and then re-expands with antiparticle 
excitations. The exponent z as usual diverges as the Fermi surface enters the oscillatory region; it 
finds its minimum at — )• 1, where it passes below z = 2. The excitation widths over energies 
also diverge before the maximum /x/j value is reached, and the corresponding pole moves off the 
physical sheet on the complex u plane. 

Case 4: T^(i'~5'i) and ^^(i'f'^l), dual to Tr A3Z3 and Tr A4Z2. 

These modes have (91,^2) = {pi,P2) = {\^-\) and (mi,m2) = {\-,-\)\ kp is given in 

figure [s] and Vkp, z and F/w in figure [oj 

The behavior in this fermion is similar to case 2, with /cqsc vanishing at ^ir — )• 0; in fact the 
charges and mass values conspire so that the first two terms in the //r — )• expansion vanish, 
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Figure 10: The values of kF/fJ'2 for case 5, Tr X3Z1 and Tr A4Z1. 
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Figure 11: The values of Vkp, z and T/u for case 5, Tr X3Z1 and Tr A4Z1. 



leaving a cuspy 0{fi^) deviation away from the limiting value kosc/ 1^-2 — ^ — 1/\/2- 

The total charge is (73 = 3/2, and again there is a single Fermi surface. For this example 
the Fermi surface never enters into the oscillatory region that we can discern, although it tracks 
it rather closely. The Fermi momentum again goes through zero, indicating a Fermi surface of 
particles transitioning to one of antiparticles. The exponent z again reaches a minimum around 
z ~ 2, and the ratio T /uj has a maximum around T /uj < 0.23, going to zero as ^ir — t- 1. 

Case 5: t^^I'"!'!) and T^(f'5'~|), dual to Tr A3Z1 and Tr A4Z1. 

These modes have (51,92) = (§,0), (pi,P2) = (-1,0) and (mi, 711,2) = (-\-,\)', kp is given in 



figure 10 and u^^, z and F/w in figure 
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In this case the oscillatory region vanishes for /i/j < 0.53 as the effective mass overwhelms the 
effective electric coupling. The effective charge is = 3/2, and there is again one Fermi surface; 
as with the other two = 3/2 cases, the Fermi surface is above the oscillatory region. As fin — )• 1, 
the Fermi momentum appears to meet the oscillatory region. It departs from the oscillatory region 
before fcosc ceases to exist, as in case 3, and interestingly, appears to head towards kp/ij,2 — >■ 0. 

This fermion is remarkable because it has 52 = and p2 = 0; it is unaware of the gauge field. 
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The effective electric coupling (^e)^^ given in (121) in general has the leading behavior as fiR — )• 



of {qe)lg —7- (j'l/4; with q2 = 0, the subleading behavior does not kick in until (^e)^^ — )• 9/xfj/8. The 
effective electric coupling is thus very close to zero for much of the range of ^r. 

The width ratio T/u reaches a maximum of T/u ~ 0.45 in the middle of the range of fiR, 
dropping off to zero on either side. Perhaps most interestingly, v^p steadily climbs as /j^r — )• 0, 
approaching z^^^ —7-1/2; this is the limit of the marginal Fermi liquid, where the non- Fermi liquid 
crosses over to an ordinary Fermi liquid. Correspondingly, the exponent z approaches z — ?■ 1. 

It is tempting to look for an explanation for such limiting behavior. As previously noted, the 
limit /Uij, — 7- is related to the BPS Coulomb branch solution of the 1-charge black hole, although 
the order of limits is subtle as ^2 remains nonzero. However, this fermion — which we call the 
"marginal Fermi liquid" fermion — has q2 = P2 = 0, and hence is unaware of ^2- Thus, as fipt — )• 0, 
its Dirac equation becomes a Dirac equation in the Coulomb branch background. In the next 
section, we derive these results analytically by matching onto the Coulomb branch solution. 

7 The 1-charge black hole and the marginal Fermi liquid fermion 

In this section we explore the BPS Coulomb branch solution that is the r^f — t- limit of the 1- 
charge black hole. We first describe how the fiji — )• limit of the "marginal Fermi liquid" fermion 
is controlled by the infinitesimal-fc limit of the Coulomb branch solution. We then describe the 
solutions for all our fermionic modes at finite k, which like other related modes explored in the 
literature manifest a continuous distribution over a gap. We then describe conductivities in this 
background, and argue that the gauge field is superconducting while is insulating, as well 
as considering the extension of the conductivities into the 2+1-charge backgrounds, where the 
superconducting delta function in the conductivity broadens to a Drude peak. 

7.1 Derivation of the marginal Fermi liquid limit 

-(- -) -(- - — i"! 

The fermions ^'^ ^iid x , referred to as case 5 in the previous section, possess the 

distinction of having q2 = P2 = 0. As a result, these modes are completely insensitive to the 
$2 gauge field, and as described before there is a limit where the Dirac equation in the 2+1- 
charge extremal background approaches that of the BPS Coulomb branch background, where things 
simplify. As we shall see, the existence of a Fermi surface approaching the marginal Fermi liquid 
value of = 1/2 can be determined analytically there. 

It is interesting to consider the behavior of kp/fii, instead of kp/ij,2, for solutions of this Dirac 
equation. For the other cases, /ii — )• with /i2 fixed is the straightforward realization of — )• 0, 
and kF/ij,2 and kosc/fJ'2 remain finite in that limit. For the MFL fermion, however, it does not feel 



/i2 and so allowing //2 — ^ 00 with ni fixed is a reasonable limit. As can be seen in figure 12, in this 



34 




Figure 12: The values of kp/lJ-i for case 5, the MFL fermions Tr X3Z1 and Tr A4Z1. 



normahzation we have the result 



(134) 



over the entire range of fiR] the line kp/ fii = 1/2 is shown explicitly. As we shall see, in the fiR — )• 
limit this relationship can be shown to be exact. It is possible it is exact for the entire range of hr, 
but we will not show this. 

The 2+1-charge extremal backgrounds, having eliminated Q2 to impose extremality, depend on 
the length scales rn and Qi. The limit approaching the Coulomb branch solution is rn ^ with 
Qi finite; thus there is a separation of scales, rn ^ Qi- Studying Fermi surfaces requires setting 
the right boundary conditions at the horizon. Consequently, we will first study the region "zoomed 
in" around rn', keeping rn finite means Qi — 00. We keep k finite in this limit as well, but it 
drops out of the equations. The Dirac equation for all four modes reduces to 



which has the solutions 



*±i/2(r) = (r2-r|,)±5. 



(135) 



(136) 



Near the horizon, these approach 



'^±i/2{r vh) [2rH{r - rn)] 



±1 



(137) 



corresponding to a value of Vk = 1/2. Indeed we may verify that the formula (110) for evaluated 
for this fermion in the large Qi limit gives z/^ — )• 1/2. 

Thus to pick the ingoing/regular boundary condition, we pick the positive exponent = +1/2. 
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Meanwhile, at large r these solutions become 

'^+i/2{r^^)^r + 0{r'jj/r), ^_i/2{r ^ oo) ^ ^ (138) 

Thus the regular solution at the horizon goes like r far away from the horizon. 

At values of r much greater than rn one finds the scale Qi, which we scaled away in the previous 
limit. Let us now "zoom in" on this region, taking Qi finite and ru 0. Here we approach the 
Coulomb branch solution (modulo the behavior of which this fermion is insensitive to) where 
the Dirac equation simplifies. Since we had k of order rn in the inner region, we will do the same 
in this outer region, setting 

k = kpTH , (139) 

and holding kp finite in the limit. 

To find a Fermi surface we study the ^'1+ mode (as usual ^2+ has the same solutions with 
A; — )• —k), and the solutions are 



r 2r(7-^ + Qi) / 

where A and B are constants. We have chosen these names for the constants for good reason, for 
in the large-r limit we find 

M/i+(r^oo)^ Ari/2 + i?5jl^r-3/2, (141) 



which are the proper boundary scalings for niL = 1/2 according to (100); A is the source term. 
For a Fermi surface, A will vanish when infalling boundary conditions are imposed. 

To determine what is implied by infalling boundary conditions, we match the small-r limit 
of this outer region to the large-r limit of the inner region just discussed, patching the solution 



together between the scales rn and Qi. The small-r limit of the outer region solutions (140) is 



^Mr^O)^A{^+ +B\-^ + ^^ . (142) 



To match to the Vk = +1/2 inner region solution, this must go like r; thus infalling boundary 
conditions requires 

, (k„ -1)B , , 

^ = \r.2 ■ (143) 



Thus for generic kp, there is a nonzero source term. However, for the special value 

kp = l k = rH, (144) 



36 



we find A = when we impose infalling boundary conditions, and thus there is a Fermi surface. 
Moreover at — )• 0, we have 

jjif^R^O)^^. (145) 

Thus we expect that as Qi/rn — t- oo, which corresponds to — ^ 0, there will be a Fermi 
surface in the family of extremal 2+1-charge black holes with parameters approaching 

kp rn = Y ^ ^kp-^^- (146) 



This exactly matches the results obtained previously. Figure 12 suggests that the relation kp ~ /ii/2 



continues to hold away from /u/j = 0, although kp = rn no longer will. 

7.2 Finite momentum fermion solutions for the Coulomb branch solution 

In the above analysis, the infinitesimal- /c limit of fermion fluctuations in the Coulomb branch 
background was matched on to the limit of the 2-|-l-charge extremal black holes for the marginal 
fermi liquid fermion. One may also consider the finite-fc fermion fluctuations, and in the true 
Coulomb branch background rather than the 2-|-l-charge extremal limit, which we discuss here. 

As discussed, this background is a zero-temperature, zero-chemical potential renormalization 
group flow geometry, referred to as the "n = 2" Coulomb branch flow in [50] , where the fluctuation 
equation for the transverse traceless graviton modes (which appear as a free scalar) were solved. 
The coupled system of the scalar and the trace of the graviton was considered in [59] , [60] , and then 
several other modes, including the spin-1/2 fermions that couple to the gravitini were analyzed in 
|51j . It is straightforward for us to do the same for the fermions we study here. 

In [50l [SH [60l El] fluctuation equations were solved subject to a regularity condition in the 
deep interior. The solutions were hypergeometrics characterized by a parameter a, defined in our 
notation as 

1 1 / tA^ 

n ' ^^^^^ 

where = uP' — k"^. The hypergeometric solutions led to two-point functions with a mass gap at 

Ai = |i , (148) 

and a continuous distribution above the gap, appropriate for a state of A/" = 4 Super- Yang-Mills on 
the Coulomb branch. 

We may perform the analogous analysis for the uncoupled spin-1/2 fields. Since the gauge fields 
are zero, the equation cares only about the mass parameters (mi, 1712), which for the moment we 
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leave arbitrary. We find for the solution regular at the deep interior (reverting to x using (91)), 
;^ = (l-^;)i+'"i+'"2^;"-i-BFQ+mi + a„i,l + mi + 2m2 + a„,,2 + 2a„i,t;^ , (149) 



where the variable v is given in (20), and we have defined 



a. 



mi 



For all our uncoupled fermions, we have \mi\ = i, and hence ami — ^ a as given in (147). For all 



the uncoupled fields in the 20, both {mi, 1712) = (1/2,-1/4) and (mi, 771,2) = (~l/2,3/4), we find 

;^20 = (1 - vy/\''-^/^F (1 + a, 2 + a, 2 + 2a, v) , (151) 

while for the fermions in the 4, with (mi, m2) = (1/2, 1/4), we have 

XA = {1- vf'\''-^''^F (2 + a, 2 + a, 2 + 2a, v) , (152) 

The X20 result exactly matches the positive-chirality solution found in the sector mixed with the 



gravitino in equation (188) of |51j : the X4 solution has the same mass gap (148) and continuous 
spectrum as well. This completes the analysis of [51j to the entire fermionic spectrum. 

We note that for a general ttii, the solution would manifest a different gap than found for all 
the other modes in [50l EHl IMl EI] • The universal spectrum is presumably a consequence of the 
underlying supersymmetry, which constrains the possible values of 777i. 

7.3 Conductivities in extremal backgrounds 

Conductivity studies complement calculations of holographic Fermi surfaces, as both should be 
interpretable in terms of A/" = 4 super- Yang-Mills theory. Here we explain how to formulate bulk 
differential equations which determine the conductivities. We will give explicit closed form results 



for the Coulomb branch solution, defined through the limit (21). The nearby extremal 2+1-charge 
black holes are the ones which exhibit holographic Fermi surfaces which are close to marginal Fermi 
liquids, and we will provide some numerical results on their conductivities as well. 

In order to study conductivities at zero wave-number but non-zero frequency, we may turn on 
gauge fields at linear order: 

a^ = e-''^%{r) A, = e-"^'B,{r) . (153) 
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In general it is necessary at the same time to perturb the metric at hnear order: 



(154) 



Other components of the gauge field may be set to zero as a gauge choice, and other perturbations 



may be decoupled at linear order from the ones indicated in (153) and (154). The perturbed 
Einstein equations reduce to 



(155) 



where X = e The equation (155) may be used to eliminate htx from the perturbed Maxwell 

equations, with the results 



h X 



e^^hX^ 



B':. + i2A' -B' + - —]B'^ ^T^6x + TTT^ Bx = 0. 



(156) 



h X J""" e^^h ^' e2^/i2 
An immediate consequence of these equations is that the flux 



^ = ie^^-^h (xXKb. + ^^B*Xb, 



(157) 



is conserved, in the sense that drJ-^ = 0. It is straightforward to show that is the imaginary part 
of the on-shell action, up to an uninteresting prefactor related to the gravitational constant. Thus 
it can be used to compute the real part of the conductivities without extracting the full two-point 
functions of the currents dual to and a^. In particular, if we define 



Aa = lim {bx,Bx) 



(158) 



then the Green's functions of the dual currents jx and Jx satisfies 



J" - ^* (Im G^^ ,^p)Ai3 , 



(159) 



where by ImG'^ we mean ^{G^ — G^^), and in defining G^ we insist upon purely infalling boundary 
conditions at the horizon. The matrix of conductivities satisfies 



Re ansii^) = Im 



(160) 



where by Re a we mean ^{a + ci^). Because Reo" is hermitian, it is diagonalizable. Positivity 
properties of the Green's function guarantee that these eigenvalues are positive. The imaginary 
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part of the conductivity, corresponding to the real part of G^^, has an additive ambiguity, Hnear 
in uj, which can only be resolved through a suitable holographic renormalization prescription. Our 
approach of side-stepping this difficulty by using a conserved flux to compute only the imaginary 
part of the Green's function is similar to the one in [61]; however, we will see that some useful 
information about the real part of the Green's function is available without the full technology of 
holographic renormalization. 

7.4 Conductivities for the Coulomb branch solution 



In the Coulomb branch limit defined in (21 ), the equations of motion for htx, bx, and Bx all decouple 



from one another. The bx and Bx equations can be simplified to 

3 oo^L^ 

In fact one equation can be transformed into the other by the rescaling 

bx = X-^Bx , (162) 

so the solutions to one may be derived from the other; the general solutions to these equations may 
be expressed in terms of hypergeometric functions. These equations were studied in [5T] , where it 
was shown that the solutions are related to the solution of the massless scalar fluctuation; however, 
it is useful to review the relevant parts of the calculation from a perspective that will facilitate the 
numerical work we will do in the next section. In the rest of this section will be convenient to set 



L = Qi = 1, so also for the gap (148) Ai = 1. 



In order to determine which solutions are physical, we must examine the behavior near r = 0, 



where the independent solutions to (161) take the form 

for r->-0, (163) 
where we have introduced an ie prescription so that the desired solution corresponds to choosing 



the plus sign in each exponent, associated to the less singular solution. The full solutions to (161) 
are then 

bx = jg / (1 + r (2 + a.) ^2+2a. (1 ^ 2 + a^- 2 + 2a.; -r') 
r (2 + 2a.) 



B. = i32 ^^^^"rl^^J^t""^ r2»" 2F1 (a., 1 + a.; 2 + 2a.; -r^) 
r (2 + 2a.) 



(164) 
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where the Ba are integration constants, and 

a^^a{k = 0) = -^ + ^^/l^. (165) 
The gamma functions could have been soaked into the definitions of Ba, but including them ex- 



plicitly makes the expansion near the boundary simpler. In (164) and below, we mean to continue 
with the prescription oj ^ uj + ie, but for simplicity of notation we do not write it explicitly. Noting 
that the hypergeometric functions approach unity as r — )• 0, it is easy to see that 

J^={\b,\' + 2\B,\^)e{io^-l) Va;2-1. (166) 

Using standard identities for hypergeometric functions, one may show that 

Aa = MapBp where M^p = ^ , tw9 ^ ^ il !l ' (^^T) 

The matrix Map is convenient because it allows us to quickly extract the real part of the conduc- 
tivities from the flux: 

1 d'^F 

Rea„^(^) = -(M-it)^^__(M"i)5^. (168) 



Because the equations of motion (161) decouple, Map and therefore ^euapioj) are diagonal. The 



conductivity eigenvalues are found to be 

Rec7a(w) = ^Rea^(^) = ^0(^2_-^)^^^l^W^Hl for / . (169) 

Again we see the structure of a continuum above a gap; this is characteristic of both insulators 
and hard- gapped superconductors. The difference between the two, from the point of view of 
conductivities, is that superconductors have an additional contribution to Re(Ta(a;) proportional to 
5{ijj) (infinite DC conductivity), whereas insulators do not. 

We expect that the Coulomb branch state is an insulator with respect to and a supercon- 
ductor with respect to a^. The reason is that, as demonstrated in [53] (and further discussed in 
|50|). the Coulomb branch state under consideration consists of a uniform disk of D3-branes, spread 
out in the Xi and X2 directions, but localized to the plane = X4 = X5 = Xq = 0, where Xi 
are the six real scalars of A/" = 4 super- Yang-Mills. Rotations in the X1-X2 directions correspond 
to the U{l)i gauge symmetry of in supergravity, while simultaneous rotations in the X3-X4 and 
X^-Xq directions correspond to U{1)2 and A^j. U{1)2 is genuinely a symmetry of the Coulomb 
branch state because each of the D3-branes lies at a fixed point of the corresponding rotations in 
ten dimensions. U{l)i is also a symmetry, but only insofar as we regard the disk of branes as truly 
uniform. At finite A^, assuming individual D3-branes must take definite positions within the disk. 
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U{l)i is explicitly broken. Even in the formal limit of infinite A^, corresponding to calculations 
in classical supergravity, it is reasonable to expect that the main features of an explicitly broken 
U{l)i symmetry will arise, namely infinite DC conductivity and the Meissner effect. 

The most efficient way to check the expectations explained in the previous paragraph is to 



use the wave-functions (164) to extract enough of the Green's function to see the infinite DC 



conductivity. To this end, consider the expansions 



S.X = Bo 



(170) 



Formally, we may compute the Green's functions through the expression 



lim T'^ = AlG^^^^pAp 



(171) 



where 



„2A-B 



(172) 



The expression (171) essentially follows from considerations of the on-shell action, i.e. the standard 



AdS/CFT Green's function prescription [5l|6]; see [62l!6l]. But it is a formal expression due to the 



in the expansions (170), which gives rise to a logarithmic divergence proportional to cj^ in the 



limit on the left hand side of (171). In holographic renormalization |63], this divergence is canceled 



through the introduction of a boundary counterterm proportional to +2F^„. The counterterm 
only effects the real part of the Green's function, and only at quadratic order in uj. Thus we do not 
need to establish a definite renormalization prescription if our aim is to read off the w-independent 



term. In short, by plugging (170) into (171) and ignoring the divergence, we obtain 



G'^x,a = -2 + O(w^) G^^^A = C>{i^'^) at smah cj and = 0. 



(173) 



Recalling the familiar relation. 



a{uj « 0) = = G^(0) i-iV- - t:5{u) ) , 



(174) 



and putting (173) together with (169), we arrive at 



Re aaioj) = 2tt5{uj) + ^^(^^ - 1) tanh — 



Re(JA(w) = t:oj6{uP' — l)tanh 



vrVw^ — 1 



(175) 
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confirming tliat tlie Coulomb branch solution has infinite DC conductivity for the current dual to 
and zero DC conductivity for the current dual to A^. 

To understand the Meissner effect, we first recall that the two-point function for a conserved 
current must take the form 

Gmnik) = Gik') (^n^n " ^) , (176) 

where k'^ = g^^kmkn and we choose mostly plus signature for the boundary metric Qmn-, which we 



assume to be the flat Minkowski metric. Putting (173) and (176) together, we find 



Gf(0) = -2 G^{Q) = Q. (177) 

We can now check that there is a Meissner effect for and not A^, using the London equation. 
Following arguments of [61] (see also the discussion in [65]), one finds at long wavelengths that 
j(0, k) = G^(0)a(0, k), or in position space j{x) = G^{0)a{x), provided we make the gauge choice 
V • a = 0. Assuming also the Maxwell equation V x Ba = j for the magnetic field i3„ = V x a, we 
arrive at 

[V2 + Gf(0)]a = 0. (178) 

Because G^(0) < 0, solutions in the presence of the condensate are exponentially rising in some 
direction and falling in another. The resolution is that the magnetic field penetrates into the 



condensate a characteristic distance l/y^— G^(0) and destroys the condensate where it is large. 
It may seem there is a loophole in the argument of the preceding paragraph: in computing 



G^{0) and G^(0), we first set /c = and then took w 0. But in (jlTSj) we used G^(0) in a 
different limit: uj = first with k made small afterward. However, due to the four-dimensional 
Lorentz invariance of the system, fluctuation solutions with nonzero k as well as co just take the 



form (164) with a^i — )• a, depending only on the combination = uP' — . The limit is then smooth 



regardless of whether w or /c is taken to zero first. 

7.5 Conductivities in the marginal Fermi liquid regime 

We now wish to investigate conductivities of extremal 2+1-dimensional black holes in the regime 
where marginal Fermi liquid behavior was found, namely large Qi/th, corresponding to ^.r — t- 0; 
we may think of this as "doping" the 1-charge Coulomb branch background with a little Q2- The 



equations of motion (156) are complicated enough that it would be unenlightening to reproduce 
them here; however recall the key feature that mixing between and yl^ is proportional to ^'i^^, 
and so it vanishes in the strict r// — )• limit but is non-zero at any finite r^. For our class of 
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extremal black holes, the coupled equations (156) may be diagonalized by either of the two choices 

V2- 



(179) 



where it is clear that near fiji — )• these correspond to the "mostly B^^ and "mostly ft^^" solutions, 
respectively. 

Another key difference between the Coulomb branch case and the case of extremal black holes 
with small but non-zero is that there is an AdS2 x region close to the horizon. The radial 
dependence of the perturbations near the horizon may be expressed as 



i(l + 3r'jj+4rjj)ui 



(180) 



B^ = ^2e^(^^^^^(^+^ (r - rn) ^V^^'+^^%^' [1 + 0{r - rn)] ■ 



The leading behavior is the same (up to the integration constants Ba), but the corrections in square 



brackets are not. It is worth noting that the asymptotic forms (180) do not reduce to solutions of 
the form (163) when rn — )• 0. This is because of order of limits: In (163) we set rn = first and 



then took r — )• 0; here we take r close to the horizon at finite rjj first. 

Conceptually, our strategy is very similar to the one outlined in the previous subsection: We 



first use (180) to write the flux as 



UJ 



2r2 



(181) 



and then we relate the coefficients Ba to the boundary values Aa of the gauge fields in preparation 
for computing the Green's function. Our numerical scheme is to start with the asymptotic forms 



(180), including explicit expressions for the O^r — rn) corrections, and use them to seed a numerical 
integration of the equations of motion out to large r, where we evaluate Aa- 

The linearity of the equations guarantees that Aa = MajsBp for some M^^, but because the 
equations couple, Man is not diagonal. It can however be extracted from numerical solutions for 



any given uj. Then we extract the matrix of conductivities precisely as indicated in (168), using the 



near-horizon form (181) for the flux. Finally, we extract conductivity eigenvalues, corresponding 



to the diagonalization (179). The results are shown in figure 13 



Several points are worth noting: 

• As r// — )• 0, the mixing between a^. and becomes small , and as a result, the eigenvectors 
of Reo"Q,^(cj) are nearly | j and f j, consistent with (|179|). We will refer to the corre- 



sponding conductivity eigenvalues as (Ta{oj) and (Ta(<^), since they predominantly refer to the 
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Figure 13: The eigenvalues of the real part of the conductivity matrix for the currents dual to Ux 
and Ax, as functions of w, for black hole backgrounds close to the Coulomb branch solution. Red 
dots show evaluations of the conductivities for rn = 1/10; blue dots are for rn = 1/100; and the 



green curves are the analytic expressions (169) for the Coulomb branch solutions. The blue and 



red curves are based on the Drude model, as explained in the main text. 



conductivities of the gauge fields dual to an and Au, respectively. 



• For a; > 1, and cryi(cj) converge quickly to the results (169) found for the Coulomb 
branch. In particular, (Ja{uj) |cTyi(w). 

• For Lu < 1, aA ^ o"a- Thus aA and aa cross at w f« 1. 

• (TA ~ io'^rfj for small u and small r^- For fixed rn, this power law behavior is approximately 
respected all the way up to w 0.3. 

• For small cu and small rn, one finds 



1 + wV^ " 8rj^ 



Again the fit is good up to a; 0.3. 



The functional form in ( 182 ) is based on the Drude model, in which more generally 

ReCJDrude = ."no , (183) 

and, conventionally, r is the scattering time while do = ne^r /m is a constant related to the density 



of charge carriers. We note that for small fio, the form (182) converges in the sense of distributions 
to a a ~ TT6{ijj) at small w, which is half as large as the Coulomb branch result ( 175| )p| This apparent 
®In taking this limit, it is important that we are considering uj to be valued over the entire real hne. Alternatively, 
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discrepancy is surprising, but it may find its explanation in the existence of an additional Tr6{uj) 
contribution to cTa(w) which our present methods do not allow us to detect. 



8 The 2-charge black hole and soft gap 

We turn now to a study of fermion fluctuations in the 2-charge black hole, which has similarities 
and differences to the 1-charge case; a difference is that here, there is a genuine extremal black hole 



with a charge density, breaking four-dimensional Lorentz invariance. As described in section 2.4 
like the 1-charge case the limit of the extremal 2+1-charge black holes is not precisely the 2- 
charge extremal black hole; unlike the 1-charge case, however, we have no fermion fluctuations that 
ignore the discrepancy. Therefore the behavior of our fermion fluctuations in the extremal 2-charge 
solution is disconnected from the 2+1-charge results of section |6j We can nonetheless study them 
in their own right. The near-horizon region of the 2-charge black hole is not AdS2 x M^, but is 
conformal to it l66j and the ten-dimensional lift contains an AdS^ factor 1201 . 



8.1 Fermion eigenvectors in 2-charge background 

The 2-charge black hole has = 0, and hence the charges qi and pi are irrelevant. A number of 
Dirac equations that are distinct for the 2+1-charge background now coincide. Below we list the 
positively-charged modes, both from "maximal" eigenvectors and the gravitino-free non-maximals: 





mi 
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1 
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The first column counts the number of modes of each chirality with the corresponding set of charges. 
We note that each charge is chirally balanced by itself, unlike in the 2+1-charge case; this reflects 
the fact that U{1)2 is non-anomalous. 

There are also negatively-charged conjugate modes with opposite q2 and p2- Finally, there are 
neutral modes, with balanced chiralities and two different mass functions: 





nil 


7712 


Q2 


P2 


2x, 


2x 


1 

2 


3 

4 








2x, 


2x 


1 
2 


1 

4 









if we restrict lj > 0, then the distribution 5{ijj) should be regarded as supported half in the cj > domain and half in 
the a; < domain. In this alternative approach we would say dw n5{uj) = ^, and also limCT„_+o duj -^^'^°2^-2 = f 
where e stays fixed in the limit. 
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8.2 Near-horizon analysis: extremal 2-charge case 



We can perform a general near- horizon analysis of the Dirac equation ( 89 ) in the extremal 2-charge 
case. Near the horizon we have 



.2A 



.25 



L2 



while for the scalar, 



6 \ r 



^2/3 ^2 



2r^ 
01 



$9 



r 

'2Q2L 



(184) 



(185) 



Like the 1-charge case but unlike the (2+l)-charge cases, the scalar diverges at the horizon. The 
mass term in the limit becomes 

^(^) ^ ?!^^ , (186) 



Lr4/3 ' 

which is also divergent. Looking at the second order equations, we find 



4^ 



a± 



0, 



(187) 



The two dominant contributions to the zero-derivative term are the frequency and the diverging 
mass. For a; 7^ this has the solutions 



^ ~ exp 



(188) 



for either ^^0+ or We notice that for infinitesimal cj, these solutions are not propagating 

waves; the divergence of the scalar has created a certain minimal value of necessary to make an 
infalling solution: 



A2 



2V2|m2|Q2 
L2 



(189) 



We can also study the u; = equation directly, and then take the small r limit, The different order 
of limits changes the coefficient of the first derivative term: 



2/o2 



r 



2m^Q2 



0. 



(190) 



Notice how this equation differs from the analogous case for the 2+1-charge black holes (109) in 



having a leading 1/r^ term generated by the diverging mass, dominating over the l/r^ term that 



would contain Equation (190) has modified Bessel function solutions 



:i9i) 



47 



which for small r have the behaviors 



r^/^ exp 



(192) 



similar to the finite-w solutions ( 188 ) with an additional r-*^/-^ factor 



The behavior brought on by the diverging mass in these known J\f = 8 gauged supergravity 
modes should be compared to [66], where Dirac equations with constant mass were studied in the 



extremal 2-charge black hole; there no analog to (189) was found. This is a case in which the 



top-down example is instructive, as the simpler solutions of [66] are apparently not present in any 
of the actual modes known to be dual to operators = 4 Super- Yang-Mills; it would be interesting 
to diagonalize the gravitino sector to see if the behavior in the solutions of [66j is present there. 

8.3 Fermi surfaces for 2-charge black hole 



We numerically solved the Dirac equations listed in section 8.1 at uj = 0, subject to the regular 



boundary condition at the horizon, which is the Kq Bessel function in (191), leading to the minus 



sign in (192). As in the 2+1-charge case, we solved the equations out to the boundary and looked 



for values of k = kp such that A{kF) vanishes. 

For the first two equations in the table, both with mass values (mi, m2) = (1/2, —1/4), we found 
such Fermi-surface-type singularities. In the first case, that with {q2,P2) = (2,0) the corresponding 
value of kp is 



kp 
^J'2 



0.83934. 



while for the second case of {q2,P2) = (1, —1/2) we obtain 

kp 



0.05202. 



M2 



(193) 



(194) 



We find no Fermi surface singularities for the remaining equations. Unlike the (2+l)-charge case. 



there is no analogue of uj^ and no analogous excitations; the form of the fluctuations ( 188 ) prevents 
modes with small uj and k± = k — kp from manifesting an ordinary infalling mode. To find such a 



solution one must go to uj > A2, where the small-o; form of the Green's function (128) no longer 



obtains. A more thorough understanding of these fluctuations is naturally of interest, and we leave 
it to future work. 



8.4 Conductivities for the 2-charge black hole 

One may also consider the gauge field fluctuations and associated conductivities. As in the 1-charge 



background, the bx and Bx fluctuation equations (156) decouple, and one can be transformed into 
the other by (162). While we cannot solve the equations analytically for general uj, the a; = 
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Figure 14: The real part of the conductivity for the currents dual to bx and B^, as functions of oj, 
for the extremal 2-charge black hole. The blue (solid) curve corresponds to the mode, and has 
a delta function contribution at zero frequency. The green (dot-dashed) line is for the current dual 
to bx, and the DC conductivity is insulated. Note the presence of a soft gap in each. 



solutions are easy to obtain (the Bx case was previously found in [67] ) . The zero frequency equations 
of motion are 

^„ , 10 + 9r2 + 3r4 



b':. + {- + J^]b:. 



(r + r3)(2 + r2) °' 

r ' 2 + (2 + r2)(r + r3)2 ' 

where we have set Q2 = L = 1, and the solutions regular at the horizon look like 

^2 



(195) 



hx = Bi, Bx = B2 



(196) 



In the spirit of (172), we can formally define the complex fluxes from (161), in this case 



-h e 



2A-B / x^8 



XHldrbx + —,BldrBx 



(197) 



Inserting the regular solutions (196) into ( 171 ) via the flux ( 197), we thus extract the zero frequency 



contributions to the Green's function 



Gi(0) 



(198) 



Accordingly, the real part of the conductivity for the current dual to Bx has a delta function con- 
tribution 4:Tt6{u}), while the current dual to bx does not. Thus in both the 1-charge and 2-charge 
cases, it is the gauge field turned on in that class of backgrounds whose dual current manifests 
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infinite conductivity at zero frequency. Here, however, we have a nonzero charge density (15), and 
lack the four-dimensional Lorentz invariance that we used in the 1-charge case to demonstrate a 
Meissner effect. Thus despite the complementary appearance of the two examples, a more care- 
ful investigation is necessary to find out whether the delta function is entirely a consequence of 
translation invariance [68], or if there really is superconductivityj^ 

While we cannot analytically solve the fluctuation equations ( |195 ) for nonzero w, we are able 
to do so numerically. In figure 14, the real part of the optical conductivity for both currents is 
plotted as a function of frequency. The most notable feature in this picture is the appearance of 
a "soft" gap at Aq- = 2Q2/LP'. This is to be contrasted with the "hard" gap characteristic of the 



1-charge Coulomb branch conductivities (169). Interestingly, this gap appears to be larger than 



the scale (189) felt by the fluctuations about the 2-charge black hole Fermi surfaces by a factor 
A,/A2 = 2^/2. 



9 Conclusions 

We have shown that Fermi surface behavior is ubiquitous in strongly coupled A/" = 4 Super- Yang- 
Mills theory at zero temperature and finite density; indeed for every realized value of the chemical 
potentials there are multiple fermionic modes manifesting a Fermi surface singularity. These Fermi 
surfaces are non- Fermi liquids for generic backgrounds; indeed they are often in the deep non- 
Fermi regime, with excitations characterized by large dispersion relation exponents, as well as 
small width/energy ratios. 

There are exceptions to this general behavior. In particular, one fermion approaches the 
marginal Fermi liquid limit, as well as vanishing zero-temperature entropy. This limit is con- 
trolled by the 1-charge Coulomb branch solution; both this and the 2-charge extremal black hole 
(which also has zero entropy at zero temperature) have the conductivity associated to the active 
gauge field infinite at zero frequency, while the other is insulating; the gaps are hard and soft, 
respectively. 

There are a number of open questions, many related to the dual field theory, which thanks to the 
top-down construction one may explore explicitly. One issue is the nature of the quasiparticle-like 
excitations forming the Fermi surfaces; in some top-down models arguments that have been made 
for a "mesino" — a gauge singlet scalar /fermion bound state (see for example |32l [55]). In 
|49) . it was pointed out that the manifest A^^ scaling of the top-down Green's functions suggests a 
particle in the adjoint, coinciding with a suggestion [20] that the gaugino itself generated the Fermi 
surface. Here we see, however, that operators containing the same gaugino but different scalars 
(such as e.g. Tr A1Z2 and Tr XiZi, case 1 and case 2 respectively) do not have the same Fermi 
surface singularities at values of /i/j away from the 3-charge black hole point. Thus the scalars are 
^We thank C. Herzog for a discussion on this point. 
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clearly important; along with the N"^ scaling this suggests that adjoint scalar /fermion bound states 
may be relevant. (Operators like Tr F+Ai are never found to have Fermi surfaces.) A Luttinger 
count of the charge density generalizing [39] (see also [301 [32l [381 Hll [691 [70] ) could be helpful. 

Another interesting question is the role of thermodynamic and superconducting instabilities. 
A number of such instabilities are known to be present at nonzero temperature in the family of 
geometries we discuss (see for example [20 1 1711 172] ). It would be valuable to understand better the 
effects of these instabilities, and whether they vanish in, for example, the marginal Fermi liquid 
limit. 

These calculations can be generalized in a number of obvious ways. All three chemical potentials 
of the 5*0(6) R-symmetry could be allowed to vary independently; we expect such a generic black 
hole will be similar to our 2+1-charge cases, as the 3-charge case is. We have also not diagonalized 
the gravitino mixing; in the 1-charge Coulomb branch case the gravitino sector led to very 
similar results as the uncoupled fermions, but it is possible in the general case new phenomena may 
lurk there. One could also extend this analysis to the three-dimensional ABJM case, where one 
mode was already studied in [39], or to the six-dimensional (2,0) theory. 

Because we find a marginal Fermi liquid and small zero-temperature entropy in the limit where 
we approach the one-charge Coulomb branch solution, it is natural to inquire how closely our 
construction might come to what is needed to describe important strongly coupled condensed 
matter systems, such as cuprates near optimal doping and heavy fermion compounds. We have 
a top-down construction, so we are finally in a position to formulate questions about a known 
strongly coupled field theory (namely = 4 super- Yang-Mills) whose answers, supergravity tells 
us, include the existence of a marginal Fermi liquid. An important feature of our Coulomb branch 
construction, however, is the existence of a condensate of the adjoint scalar fields in = 4 super- 
Yang-Mills, and these are fundamental scalars, not fermion composites. The condensate breaks 
SU{N) down to C/(l)^, and the 0{N) modes associated with the unbroken U{\)^ are massless 
because they correspond to motions of individual D3-branes in the directions transverse to the 
world-volume — which is to say, directions in which the presence of the D3-branes spontaneously 
breaks translational symmetries. One can also argue that these modes are massless simply because 
of the unbroken gauge invariance associated with them; at any rate, first-principles field-theory 
arguments show that there are 0{N) massless modes. Remarkably, we learn from supergravity 
that the 0{N'^) non-abelian degrees of freedom are all gapped in the strict Coulomb branch limit 
when Qyj^jN is large: the evidence for this is the gap in the spectral measure observed, for example. 



in (175). As we depart slightly from the Coulomb branch limit, we observe the marginal Fermi 
liquid behavior at low energies and small kp, as compared to the energy and momentum scale 
set by the gap. It is reasonable to suppose that the A/" = 4 super- Yang- Mills dynamics above the 
gap can be integrated out, leading to some effective theory that controls the marginal Fermi liquid 
behavior that we see. It would be very instructive to develop field theory arguments that support 
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this picture without exphcit reference to supergravity. It may happen that at subleading order in 
N, the 0{N) massless modes give rise to interactions that smooth out the Fermi surface singularity. 
Even if the specific construction we have laid out is not directly related to specific condensed matter 
systems, a first-principles understanding of how strongly coupled field theory dynamics leads (in 
some appropriate limit) to a marginal Fermi liquid would clearly be an important advance. 
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